Transport properties and structures of vortex matter in layered superconductors 
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In this paper we analyze the structure, phase transitions and some transport properties of the vortex 
system when the external magnetic field lies parallel to the planes in layered superconductors. We 
show that experimental results for resistivity are qualitatively consistent with numerical simulations 
that describe the melting of a commensurate rotated lattice. However for some magnetic fields, the 
structure factor indicates the occurrence of smectic peaks at an intermediate temperature regime. 



I. INTRODUCTION 

The discovery of high T c superconductivity renewed 
the interest for the thermodynamic, structural and dy- 
namical properties of vortex matter. Due to the large 
temperatures available for the vortex system and an im- 
portant number of relevant parameters, as anisotropy, 
disorder and the magnitude and direction of the exter- 
nal magnetic field, the physics of these systems is very 
rich. First and second order phase transitions between a 
low temperature solid phase an a high temperature liquid 
phase have been predicted and observed experimentally. 
The low temperature phase can be either a crystalline or 
a glassy phase depending on the aiapunt and nature of 
the disorder present in the sample.Era 

All the cuprate high T c materials have in common 
a crystalline structure based on the existence of CuO 
planes. This makes all of them anisotropic materials with 
a layered structure. In the usual convention, the c-axis 
points in the direction perpendicular to the planes and 
the a and b axis are the in-plane crystalline axis. 

Concerning the study of the vortex properties in these 
materials, most of the work was devoted to the study of 
the configuration in which the external field points along 
the c-axis, particularly in the case of experimental work. 
There is however an increasing interest in the .properties 
of the system for fields parallel to the planesJ3i3 Theoret- 
ically, in many cases, the system is simulated as a highly 
anisotropic but homogeneous system. 113 The planes how- 
ever act as a strong potential for the vortices that tend 
to localize them between planes where superconductivity 
is weak. 

In a clean system and at low temperatures, for the ex- 
ternal field parallel to the ab— planes, the vortices form 
an anisotropy - distorted Abrikosov lattice. If the lattice 
is commensurate with the periodic potential due to the 
planes, the ground state of the system is a commensurate 
structure like the one schematically shown in Fig. (la). In 
this particular case, the vortex density S in the direction 
perpendicular to the planes is modulated with a period 
given by the lattice parameter of the vortex lattice. For 
this type of structures, the most salient theoretical pre- 
diction is the existence, at an, intermediate temperature 
interval, of a smectic phase.tffi In this picture the transi- 
tion between the crystalline state and the high temper- 



ature liquid state takes place in two steps. In the inter- 
mediate regime the system develops long (or quasi-long) 
range order in the direction perpendicular to the planes, 
i.e. the liquid develops some density oscillations along 
this direction as a precursor of the frozen state. This 
phase is known as the smectic phase. The existence of 
the smectic phase is unconfirmed and there are many rel- 
evant questions concerning its nature and stability. One 
of them is what happens if the field does not produce a 
commensurate structure like the one shown in the figure. 
One possibility is that the lattice locally retains its struc- 
ture and orientation and generates discommensuration to 
accommodate to the periodic potential generated by the 
planes. If the mismatch is not too large the discommen- 
surations are far apart and the physical properties of the 
system are not very affected. In this case it could be pos- 
sible to detect the existence of the smectic phase without 
tuning the external field to its commensurate value. The 
other alternative is that as the field is changed, the vor- 
tex lattice rotates and distorts - due to the anisotropic 
properties of the system- to form a new structure that 
is commensurate with the periodic potential like the one 
shown in Fig. (lb). In this paper we will show that at low 
fields, as the field changes, between two consecutive com- 
mensurate structures of the type of Fig. (la), there are a 
number of commensurate phases of the type illustrated 
in Fig.(lb). 

For the particular case of the rotated structure shown 
in the figure, the vortex density S of the lattice is the 
same between any two consecutive planes. In this case, 
as the temperature increases, the lattice could melt going 
directly from the solid to the liquid state. If this is the 
case, it would be very difficult to observe the smectic 
phase for an arbitrary chosen external field. 

The transport properties of the vortex system for the 
configuration of interest were recently measuredliil and 
the data were compared with the theory of the smectic 
phase. Experiments and theory are only partially consis- 
tent. 

In the rest of the paper we discuss, in terms of the 
London theory with a periodic potential, the stability 
of the different commensurate phases. We also present 
some numerical simulations for the resistivity in layered 
structures and analyze the corresponding low tempera- 
ture vortex configurations. We will show that the exper- 
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imental results for the resistivity are qualitatively con- 
sistent with numerical simulations. For some fields, the 
structure factor indicates the occurrence of smectic peaks 
at an intermediate temperature regime. 



(a) 



(b) 



FIG. 1. Commensurate structures at low temperatures. In 
(a) the magnetic field applied in the y-direction generates a 
non-rotated vortex lattice for which the vortex density S is 
modulated in the direction perpendicular to the ab-planes. 
In (b) the vortex lattice is rotated an angle 8, and S is the 
same for all the planes. 



II. COMMENSURATE STATES IN THE LONDON 
APPROXIMATION 

In order to analyze the low temperature structures of 
the vortex lattice we resort to the London approach in 
an anisotropic material and in the presence of a uniaxial 
potential representing the effect of the planes. Following 
the work by Campbell, Doria and Kogar£3 (CDK), the 
free energy per unit length in the direction of vortices for 
an anisotropic material is given by: 



Fq = I (h 2 + A 2 mjj curlih curljh) 



dxdz 



(1) 



where h(a;, z) is the local magnetic field in the plane per- 
pendicular to the vortices, A 2 is proportional to the av- 
erage mass M av = (M1M2M3) 1 / 3 with M k being the 
principal values of the mass tensor My and my = 
Mij/M av is the effective-mass tensor. We consider the 
external field in the y— direction (parallel to the planes) 
that coincides with one of the principal axis of the crys- 
tal, h = (0, h y , 0). Then we take m xx — m yy = mi and 
m zz = m 3 . 

For a vortex lattice, h y (x, z) is a periodic function with 
nonzero Fourier components h y (G), where G are the re- 
ciprocal lattice vectors of the vortex lattice. The free en- 
ergy is minimized with respect to h and following CDK, 
Fq is given by: 



F 
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G 



\ 2 {miG x + m 3 G z 



(2) 



where B = fori is the magnetic induction, fo is the 
flux quantum and n is the number density of vortices. 
In this expression, the summation is over all the re- 
ciprocal lattice vectors G of the vortex lattice. The 
mass anisotropy distorts the hexagonal lattice compress- 
ing it in the z— direction and expanding it along the 
x— direction. The lattice can rotate to form an angle 9 
with the x— axis as shown in Fig (lb). For this particular 
configuration, in which the magnetic field is parallel to 
one of the principal axis of the crystal, CDK shown that 
the F contribution to the free energy does not depend 
on 9. 

The above considerations are valid for continuous 
anisotropic superconductors. In layered superconduc- 
tors, the presence of planes introduces a periodic poten- 
tial that will partially break the degeneracy in 9. In a 
first order approximation, we describe the effect of the 
planes by including a periodic potential that tends to 
localize the vortices in the interplane spacing: 



(3) 



^ = ^0+ / h y (x,z)V(z)dxdz 



where V{z) is a periodic potential with the periodicity of 
the c— axis lattice parameter. The condition for commen- 
surability of the vortex lattice with the periodic potential, 
corresponding to the situation in which all vortices are 
placed at a minimum of V(z), is given by the condition 
xj = ns where Xj is the x-component of the coordinate 
of vortex j, n is an integer and s is the c-axis lattice pa- 
rameter (the interplane distance). In the reciprocal space 
this condition gives a relation between the reciprocal lat- 
tice vectors and the vector Q of the periodic potential 
V(z): 



2?r 

Q = — z =pGi + qG 2 

s 

where p and q are integers and 

27rsin(6! + tt/3) „ 2tt7 cos(6» + tt/3) „ 
Gi — — r : 777 — x 7 77; — y 



sin7r/3 



cos tt/3 



(4) 



(5) 



2tt sin 9 2irj cos 9 

G 2 = y- 777X+— 777 y (6) 

jL svo.tt/6 L COS7T/3 



with 7 2 = yrni/rn^ the anisotropy factor and L 2 = 

(2$o)/(V3B). 

This condition can be put in the form: 



tan(0) = 



pcos(7r/3) — q 



and 



B 



psin(7r/3) 
2^o7 2 sin 2 (7r/3) 



\/3s 2 (psin(0 + tt/3) - qsm(9)) 2 



(7) 



(8) 



In Fig. (2), all angles ^ 9 ^ ir/6 satisfying the condition 
for commensurability are shown as a function of the mag- 
netic induction B. In the London approximation, since 
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the free energy is 6 independent, all commensurate states 
obtained for a given magnetic induction B are degener- 
ate. As can be seen in Fig. (2), for low fields, between 
two consecutive commensurate structures of the type of 
Fig. (la) - corresponding to = or tt/6 - there are a 
large number of commensurate phases with 8^0. As 
the field increases the number of commensurate phases 
decreases and in the limit of high magnetic fields (or high 
anisotropy) there are only two undistorted lattices corre- 
sponding to 6 = and 7r/6, and at all temperatures the 
vortex density is the same between any two consecutive 
planes. In this limit, increasing the magnetic field com- 
presses the vortex lattice in the z— direction. QwrSimple 
model is not appropriate to describe this limit ]13'EJ 




FIG. 2. Angles of rotation of the vortex lattice vs. mag- 
netic field. Each point of the diagram correponds to a com- 
mensurate state of minimum energy. Magnetic field is in 
Tesla and the parameters are those typically used to describe 
YBaCuO: s — 10 A and mass anisotropy rn\jmz = 50. 

For all angles, the high anisotropy generates structures 
where the vortex-vortex distance along the plane direc- 
tions is much larger than the distance perpendicular to 
the planes, reiving rise to vortex chains. Recently Hu 
and TachikB studied the ground state configuration of 
the vortex lattice, using Monte Carlo simulations in the 
highly anisotropic XY model for large systems. They 
found structures consisting of buckling vortex chains. We 
interpret these results as rotated lattices with discom- 
mensurations due to the mismatch of the vortex lattice 
parameter with the interplane distance. In the numerical 
simulations, transverse discretization of the space could 
be another source of discommensurations. Although in 
general, we expect the discommensuratkjtn. to form an an- 
gle of 45° with respect to the planes ,EjIL3 the boundary 
conditions in a finite sample could stabilize discommensu- 
rations parallel to the planes. The buckling vortex chains 
structures are consistent with Bitter-pattern observa- 
tions, however in the experimental case twin boundaries 
may be relevant to stabilize the observed structure Ji3a 



III. NUMERICAL SIMULATIONS 

In this section we present numerical simulations for 
the transport properties and structure factor in an 
anisotropic system described by a three-dimensional (3D) 
Josephson-junction array. The model has been exten- 
sively used and is described in Refs.[17,18] and here we 
give only a brief summary of the method. 

The equilibrium physics of this system is described by 
the Hamiltpman of a three-dimensional frustrated XY 
model:N-El 

H = -Ej costs' - tp*' - A"') + 

<i,i' >€zab— plane 



- y 



cos(^ 1 - if 1 - A 11 ) 



(9) 



<i,i' >Gc—axis 



with Ej = jTjkrs (A is the penetration length in the 
ab planes and s is the interplane distance) and A" = 

(2e/hc) f. A.dl the integral of the vector potential of the 
external magnetic field from site i to site i'. The phases 
(p l (t) are defined in the nodes of a lattice and represent 
the phase of the order parameter. The thermodynamics 
of this Hamiltonian coincides with the equilibrium prop- 
erties of the 3D Josephson-junction array. 

The dynamics of the 3D Josephson-junction array is 
contained in the time evolution of the phases <p l {t). Near- 
est neighbor nodes are coupled with Josephson junctions 
characterized by critical currents /*' and normal resis- 
tances R m . The equations describing the model are 
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1 wafo* - ip l - A 1 



2eR l 



at 



) + 
„«'(*) 



(10) 



(ii) 



Equation ( JToj ) gives the current between the nearest 
neighbor nodes i and i' with phases (p l and <p l . Here 
rf l (t) is an uncorrelated gaussian noise that incorpo- 
rates the effect of the temperature. Equation ([ll]) en- 
sures the current conservation at each node and j\ xt 
is the external current applied at node i. Equations 
( |l0| ) and ( pd| ) are numerically integrated on time using 
a Runge-Kutta method. The typical time step used is 
At = .It j (tj = 4>o /2tt RqI c ) and the number of itera- 
tions is 20, 000 < N < 100, 000. Details of the numerical 
method have been presented in previous works.E3£3 The 
mean in-plane critical currents ij are larger than the 
mean inter-plane critical currents i~ by an anisotropy 
factor 7 2 ( 7 2 = il/I^ ), and ll = I c = ^£ At the 

same time, the ratio between the in-plane resistance i?" 
and the out of plane resistance R 1 - is given by I/7 2 . We 
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also consider a small amount of disorder by taking a uni- 
form distribution of critical currents of width A defined 
as 



A = 



^Jmax _|_ Jmin^j 



(12) 



where I™ ax and /™ m are the maximum and minimum 
values of the critical current in the corresponding direc- 
tions. The disorder simulated are typically < A < .1. 
As in the previous section, we take the z— direction as 
the direction perpendicular to the planes (parallel to the 
c-axis of the crystal) and the external magnetic field 
along the y— direction. The magnetic fields simulated 
arc 1/72 < / < 1/6, where / = Ba 2 /cf)Q (a is the square 
lattice period). The values of anisotropy are 1 < -f 2 < 45 
and typical system sizes are 8 < L Xl L yi L z < 64. We cal- 
culate the resistivity in the three directions by applying a 
small probe current and evaluating the average voltage. 
For example, for the resistivity p M in the ^—direction we 
drive the system with a small current 1^ (1^ = 0.01/^) 
and measure the voltage 



h d 



¥>*)> 



(13) 



Then Plx = W 



A. Transport properties 

The results for the resistivity calculated with periodic 
boundary conditions (PBC) along the field direction and 
free boundary conditions (FBC) in the other directions 
are shown in Fig. (3). These results correspond to a sam- 
ple with / = 1/12, 7 2 = 20, L x = L y = L z = 30, A = .05 
and TV = 50, 000. We found similar results for systems 
with / = 1/48,1/24,1/6, 7 2 = 5,10,25 and the same 
sizes and disorder as the ones indicated above. The gen- 
eral behavior is qualitatively similar to the experimental 
data obtained by Grigera et al. shown in the inset for 
comparison. 




FIG. 3. Resistivities in the three directions for a system 
with FBC in x and z directions and PBC in the field direction. 
The parameters are / = 1/12, 7 2 = 20, A = .05, N = 50, 000 
and L x — L y = L z = 30. p y is the resistivity parallel to the 
field, p x is the perpendicular to the field but parallel to the 
ab-planes and p z is perpendicular both to the planes and the 
field. In the inset experimental results of Ref.[ll] are shown. 

The response of the system when the external current 
is perpendicular to the planes (and the Lorenz force is 
parallel to them) is given by p z . For this geometry, the 
pinning due to the planes is unimportant and for temper- 
atures higher than a characteristic temperature Ti the re- 
sistivity p z increases rapidly. The stable phase for T > Tj 
corresponds to a liquid phase and the rapid increase of p z 
is an indication of the high mobility of vortices parallel 
to the planes. The transition from the low temperature 
phase to the high temperature liquid phase is continu- 
ous. These results, obtained in the finite system, are not 
conclusive on whether the change of behavior observed 
at Ti is a second order phase transition or a crossover 
between two different regimes. In finite systems, thermal 
activation is observed down to low temperatures and fi- 
nite size scaling is necessary in order to characterize the 
transition. Based on previous results and the discussion 
of next section, we will refer to this behavior as contin- 
uous transition and to zero resistivity when the noise in 
the voltage is larger than its mean value, corresponding 
typically to p^ < x 10~ 5 , where p = x,y, z. 

When the external current is applied along the 
x— direction, the Lorenz force is perpendicular to the 
planes and for temperatures T > Ti the in-plane pin- 
ning dominates the vortex dynamics resulting in a small 
p x . The resistivity p x decreases strongly near a charac- 
teristic temperature T* > Ti and has a tail that extends 
down to Ti. This behavior is obtained for different val- 
ues of the anisotropy and magnetic field, it qualitatively 
reproduces the experimental observation and resembles 
the prediction of the smectic phase theory. In this the- 
ory, as the temperature is lowered, the system undergoes 
a transition to a smectic phase at a temperature T s . In 
the temperature interval Tj < T < T s the system de- 
velops long ( or quasi-long) range order in the direction 
perpendicular to the planes and the resistivity shows, in 
the vicinity of the smectic temperature T s , a critical be- 
havior of the form p x (T)- p x (T s ) cx IT-T^ 1 '" where a 
is the specific heat critical exponent. At the temperature 
Tj, long range order in the direction parallel to the planes 
appears giving rise to a crystalline structure. Since in fi- 
nite size systems the resistivity is not necessarily a good 
quantity to characterize a phase transition (in particular 
it will not show an infinite slope) we search for some sig- 
nature of a smectic phase by analyzing the structure of 
the liquid for temperatures T >Ti corresponding to the 
regime where p x shows a tail. The details of the structure 
factors calculated with PBC at different temperatures are 
presented in next section. 

Finally, we have also calculated the resistivity p y corre- 
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sponding to a transport current along the field direction 
for which there is zero Lorenz force. In this geometry, the 
resistivity goes to zero at a temperature T p > Tj. This 
behavior is due to finite size effects; dissipation occurs 
when vortices entangle to form a structure that perco- 
lates in the directions perpendicular to the external field. 
This happens at a temperature T p that is sensitive to the 
thickness (L y ) of the sample.-the thicker the sample, the 
lower is T p . Previous studieslia of the resistivity in the di- 
rection of the field show that in the thermodynamic limit 
T p coincides with Tj as experimentally observed and we 
expect the resistivity along the three directions to vanish 
at the same temperature. 
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FIG. 4. Resistivities in the three directions with PBC for 
a system with / = 1/24, j 2 = 9, L x = L z = 48 and L v = 8. 
p x , p y and p z are denned in Fig. (3). In the inset the I-V 
characteristics for the current along the three directions are 
shown. 

In order to discard boundary effects due to surface bar- 
riers, we have also calculated the resistivities and struc- 
ture factors using PBC in the three directions. Results 
for the resistivity in a highly anisotropic sample with 
weak disorder are shown in Fig. (4). For this case, the 
parameters are / = 1/24, j 2 = 9, L x — L z = 48, 
L y = 8 A = and N = 30,000. We have also done 
some simulations with PBC in the three directions for 
/ = 1/72,1/48,1/8, 7 2 = 1 - 40 and N = 10 4 - 10 5 , 
and we found similar results. While p x and p y behave 
essentially as in the previous case, as the temperature de- 
creases p z saturates at a value different from zero, which 
increases with anisotropy. The I-V characteristics for the 
current along the three directions at low temperatures are 
shown in the inset and indicates a very small value of the 
critical current in the z— direction. The saturation of the 
resistivity at low temperatures is due to the fact that the 
run was done with a value of the transport current j ex t 
larger than the critical current j c and the results corre- 
spond to a flux flow regime. A systematic study of the 
resistivities with j ex t < j c requires much more statistics 
and generates larger numerical errors. In any case, these 
results show two important points: first, for the systems 



under consideration, p x and p y are not very sensitive to 
the boundary conditions, and second, when vortices flow 
parallel to the planes, pinning is very weak. In the case 
of FBC the surface barriers generate some extra pinning 
that increases the critical current. In the thermodynamic 
limit, any rigid structure is pinned by impurities and sur- 
face effects and in this sense, it may be more appropriate 
to compare the case of FBC with experiments if results 
with j ext < j c are not available for the case of PBC. As 
a side comment, notice that the noise in p z disappears 
at a temperature of the order of Tj at which p y van- 
ishes. A simple interpretation of this effect is that in the 
flux flow regime an ordered structure generates less noise 
than a liquid. However a systematic study of the noise, 
that could give information on the nature of the different 
phases, is needed to reach definitive conclusions. 

B. Structure factor 

In this section we present results for the structure fac- 
tor calculated for different temperatures. At high tem- 
peratures, in the liquid phase, the system is highly disor- 
dered and presents vortex loop excitations. An instanta- 
neous picture of the vortices crossing a plane perpendic- 
ular to the external field (y — L y /2) is shown in Fig. (5a) 
for the same system as Fig. (4). 




(a) x (b ) 

FIG. 5. Instantaneous configuration of the vortices in a 
plane perpendicular to the field for the system of Fig. (4). 
Black dots are vortices and white ones are antivortices. In 
(a) the temperature is ^- — 1.1 and in (b) = 0. 

The vortex-antivortex pairs are small loops confined 
between two ab— planes and cut by the y = L y /2 plane. 
As the temperature decreases the system evolves towards 
an ordered solid structure. It is known that, due to nu- 
merical limitations, in three dimensions it is very difficult 
to cool the system into an ordered lattice. The type of 
structures obtained by slowly cooling the system is shown 
in Fig (5b) . We obtain this kind of structures for systems 
with / < 1/12. For higher magnetic fields, a triangular 
vortex lattice was found like the one predicted by Ref . [3] . 
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Although the obtained structure is very disordered, the 
tendency to form vortex chains (at approximately 45° in 
the figure) can be observed. This tendency, is reflected 
by the structure factor S(q±,y) defined asE3: 



5(qi,y) = 



1 



,«u. 



where N is the total number of vortices, rjj is the vortex 
charge: 1 for vortices and -1 for antivortices, and rj are 
the vortex coordinates in the direction perpendicular to 
the external field at a plane with coordinate y. We define 
£(q_i_) as the average over y of S(q±,y). 

In Fig. (6) the high and low temperature structure fac- 
tor for the system of Fig. (5) are shown. 



observed. This peak indicates the presence of a smectic 
phase. As the temperature is lowered from high temper- 
atures, first the intensity of the smectic peak increases, 
goes through a maximum and then decreases as the crys- 
talline peaks corresponding to a rotated structure in- 
crease. The wave vector g° may also be weakly tern- 

(14) 

perature dependent: The first smectic peak observed has 
qy ~ Q/3 (Q = 2-k/s) and as the temperature is lowered 
it shift towards Q/2. In Fig. (8) the temperature depen- 
dence of the smectic peak and one of the crystalline peaks 
are shown. The oscillations in the amplitude of the smec- 
tic peak are probably due to finite size effects, since only 
some discrete values of a® are consistent with the PBC. 
At temperatures T ~ Tj both the smectic and the crys- 
talline peaks coexist, indicating the coexistence of the 
two phases. 
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(b) 



FIG. 6. Structure factor of the configurations of Fig. (5). 
(a) and (b) correspond to the same temperatures showed in 
that figure. 

At high temperatures only a background of the form 
cos(q x a), where a is the lattice parameter of the junc- 
tion network, is obtained. This background is due to the 
presence of small loops which in general are confined be- 
tween to consecutive ab— planes, i.e. vortex-antivortex 
pairs which are oriented in the .T-dircction and bound at 
a distance a. At low temperatures, the background dis- 
appears and well defined peaks are observed. The peaks 
in the structure factor correspond to a rotated vortex 
lattice with a small angle 8. In Fig. (6b) we draw one 
of the corresponding rotated hexagonal reciprocal space 
unit cell. We also see that there are peaks correspond- 
ing to the same lattice structure reflected in 180° with 
respect to the z— axis. The coexistence of these two re- 
flected structures could explain the chain-like ordering 
observed in real space configurations (Fig. (5b)). For this 
particular value of the external field, we follow the evo- 
lution of the structure factor as the system is cooled and 
observe some indications of a smectic phase at interme- 
diate temperatures. A sequence of the structure factor 
at intermediate temperatures is shown in Fig. (7). 

At temperatures T > Tj, where the resistivity p x 
presents a small tail, a peak at q = (0, g°) is clearly 




FIG. 7. Structure factor for the system of Fig. (5) at in- 
termediate temperatures. A smectic peak can be observed 
at q y ~ Q/2 at T — 0.53 (a), which decreases at a lower 
temperature T=0.30 (b). 



Since the symmetry group of the smectic phase is not 
a subgroup of the one corresponding to the symmetry 
of the rotated lattice, we expect a first order transition 
between these two phases. This is consistent with the 
observation of a coexistence of the two phases in the nu- 
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merical simulations. 

For other values of the parameters studied, (/ = 1/24 
and 7 2 = 2.56,25), for which the low temperature phase 
corresponds to a rotated crystal with a larger rotation 
angle, we do not observe indications of a smectic order in 
the liquid phase. In these cases, our preliminary results 
indicate a continuous transition from the liquid to the 
frozen state without any precursor of long or quasi long 
range order in the direction perpendicular to the planes. 
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FIG. 8. Temperature dependence of the intensity of the 
smectic peak described in the text (open symbols) and a crys- 
talline peak, (q°,q°) ~ (Q/12,Q/4) (filled symbols). 

IV. SUMMARY AND DISCUSSION 

We have studied the problem of the low temperature 
structure, the thermodynamic and transport properties 
of vortices when the external field is applied parallel to 
the CuO planes in high T c superconductors. A simple 
analysis based on the London theory in the presence of a 
periodic potential indicates the possibility of a variety of 
commensurate structures that essentially correspond to 
anisotropy-distorted Abrikosov lattices rotated to match 
the periodic potential. For low fields a large number of 
degenerate structures corresponding to different rotation 
angles are obtained. Using parameters like those typi- 
cally used to describe YBaCuO, s = 10 A and a mass 
anisotropy 7711/7713 = 50, for magnetic induction between 
5 and 6 Tesla, we obtain about seven rotation angles 
which generate commensurate structures. 

The numerical simulations in a XY-model clearly show 
the tendency to form rotated structures. Our results 
shows structures with grains of twin phases. 

The transport properties qualitatively reproduce the 
experimental results. In particular, p x shows a rapid de- 
crease and a tail that extends down to T, . This behavior 
is similar to the predictions of the smectic phase theory 
although transport properties seem not to be enough to 
prove the existence of this phase at an intermediate tem- 
perature range. 



When the rotated angle of the low temperature struc- 
ture, which depends both on the value of the external 
field and the anisotropy, is small, we observe well defined 
peaks corresponding to a smectic phase. The vector Q 
and the amplitude of the smectic peak are temperature 
dependent and a first order transition is expected from 
the smectic to the crystalline phases. For the parame- 
ters that give a larger rotation angle of the crystalline 
phase, we observe a direct evolution of a system from 
the anisotropic liquid to the vortex lattice without any 
indication of an intermediate smectic phase. 
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